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I 1 Detecting heat from minor planets in the outer solar system is challenging, yet it is the 

most efficient means for constraining the albedos and sizes of Kuiper Belt Objects (KBOs) and 

■ their progeny, the Centaur objects. These physical parameters are critical, e.g., for interpreting 
spectroscopic data, deriving densities from the masses of binary systems, and predicting 
occultation tracks. Here we summarize Spitzer Space Telescope observations of 47 KBOs and 
Centaurs at wavelengths near 24 and 70 /im. We interpret the measurements using a variation 
of the Standard Thermal Model (STM) to derive the physical properties (albedo and diameter) 
of the targets. We also summarize the results of other efforts to measure the albedos and sizes 
of KBOs and Centaurs. The three or four largest KBOs appear to constitute a distinct class in 
terms of their albedos. From our Spitzer results, we find that the geometric albedo of KBOs 
and Centaurs is correlated with perihelion distance (darker objects having smaller perihelia), 
and that the albedos of KBOs (but not Centaurs) are correlated with size (larger KBOs having 
higher albedos). We also find hints that albedo may be correlated with with visible color (for 
Centaurs). Interestingly, if the color correlation is real, redder Centaurs appear to have higher 
albedos. Finally, we briefly discuss the prospects for future thermal observations of these 
primitive outer solar system objects. 



1. INTRODUCTION 

The physical properties of Kuiper Belt Objects (KBOs) 
remain poorly known nearly 15 years after the discovery 
of (15760) 1992 QBi (Jewitt and Luu, 1993). While KBOs 
can be discovered, their orbits determined, and their visible- 
light colors measured (to some extent) using modest tele- 
scopes, learning about fundamental properties such as size, 
mass, albedo, and density remains challenging. Determin- 
ing these properties for a representative sample of TNOs 
is important for several reasons. Estimating the total mass 



of material in the transneptunian region, and relating vis- 
ible magnitude frequency distributions to size- and mass- 
frequency is uncertain, at best. Quantitative interpreta- 
tion of visible and infrared spectra is impossible without 
knowledge of the albedo in those wavelength ranges. Size 
estimates, when coupled with masses determined for bi- 
nary KBOs (see Noll et al. chapter), constrain the density, 
and hence internal composition and structure, of these ob- 
jects. All of these objectives have important implications 
for physical and chemical conditions in the outer proto- 
planetary nebula, for the accretion of solid objects in the 
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outer Solar System, and for the collisional evolution of 
KBOs themselves. Of course, there is a relative wealth of 
information about Pluto and Charon, the two longest known 
KBOs, and we do not address their properties further here. 

The Centaur objects, with orbits that cross those of one 
or more of the giant planets, are thought to be the dynam- 
ical progeny of KBOs (e.g. Levison and Duncan, 1997; 
Dones et al. chapter). The Centaurs are particularly inter- 
esting both because of their direct relation to KBOs, and 
also because their orbits bring them closer to the Sun and to 
observers, where, for a given size, they are brighter at any 
wavelength than their more distant relatives. Because of 
their planet-crossing orbits, the dynamical lifetimes of Cen- 
taurs are relatively short, typically a few Myr (e.g. Horner 
etal, 2004). 

The sizes of some KBOs and Centaurs have been de- 
termined by a variety of methods. Using HST, Brown and 
Trujillo (2004) resolved the KBO 50000 Quaoar, placed an 
upper limit on the size of Sedna (Brown et al. 2004), and 
resolved 136199 Eris (Brown et al, 2006). Recently Ra- 
binowitz et al. (2005) placed constraints on the size and 
albedo of 136108 (2003 EL 6 i) based on its short rotation 
period (3.9 hr) and an analysis of the stability of a rapidly 
rotating ellipsoid. Trilling and Bernstein (2006) performed 
a similar analysis of the lightcurves of a number of small 
KBOs, obtaining constraints on their sizes and albedos. Ad- 
vances in the sensitivity of far-IR and sub-mm observato- 
ries have recently allowed the detection of thermal emis- 
sion from a sample of outer solar system objects, providing 
constraints on their sizes and albedos. Jewitt et al. (2001), 
Lellouch et al. (2002), Margot et al. (2002, 2004), Altenhoff 
et al. (2004), and Bertoldi et al. (2006) have reported 
submillimeter-millimeter observations of thermal emission 
from KBOs. Sykes et al. (1991; 1999) analyze Infrared As- 
tronomical Satellite (IRAS) thermal detections of 2060 Ch- 
iron and the Pluto-Charon system, determining their sizes 
and albedos. Far-infrared data from the Infrared Space Ob- 
servatory (ISO) were used to determine the albedos and di- 
ameters of KBOs 15789 (1993 SC), 15874 (1996 TL 66 ) 
(Thomas et al., 2000) and 2060 Chiron (Groussin et al., 
2004). Lellouch et al. (2000) studied the thermal state of 
Pluto's surface in detail using ISO. Grundy et al. (2005) 
provide a thorough review of most of the above, and in- 
clude a sample of binary KBO systems with known masses, 
to constrain the sizes and albedos of 20 KBOs. 

Spitzer Space Telescope (Spitzer hereafter) thermal ob- 
servations of KBOs and Centaurs have previously been 
reported by Stansberry et al. (2004: 29P/Schwassmann- 
Wachmann 1), Cruikshank et al. (2005: 55565 2002 AWi 97 ), 
Stansberry et al. (2006: 47171 1999 TC 36 ), Cruikshank et 
al. (2006), Grundy et al. (2007a: 65489 2003 FXi 28 ) and 
Grundy et al. (in preparation: 42355 2002 CR46. Here we 
summarize results from several Spitzer programs to mea- 
sure the thermal emission from 47 KBOs and Centaurs. 
These observations place secure constraints on the sizes 
and albedos of 42 objects, some overlapping with determi- 
nations based on other approaches mentioned above. We 



present initial conclusions regarding the relationship be- 
tween albedo and orbital and physical properties of the 
targets, and discuss future prospects for progress in this 
area. 

2. THERMAL MODELING 

Measurements of thermal emission can be used to con- 
strain the sizes, and thereby albedos, of un-resolved targets. 
Tedesco et al. (1992; 2002) used Infrared Astronomical 
Satellite (IRAS) thermal detections of asteroids to build a 
catalog of albedos and diameters. Visible observations of 
the brightness of an unresolved object are inadequate to de- 
termine its size, because that brightness is proportional to 
the product of the visible geometric albedo, py, and the 
cross-sectional area of the target. Similarly, the brightness 
in the thermal IR is proportional to the area, and is also a 
function of the temperature of the surface, which in turn 
depends on the albedo. Thus, measurements of both the 
visible and thermal brightness can be combined to solve 
for both the size of the target and its albedo. Formally the 
method requires the simultaneous solution of the following 
two equations: 
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where F is the measured flux density of the object at a 
wavelength in the visible ("vis") or thermal-infrared ("ir"); 
Fo jV i s is the visible-wavelength flux density of the Sun at 
1 AU; r and A are the object's heliocentric and geocen- 
tric distances, respectively; R is the radius of the body (as- 
sumed to be spherical); py is the geometric albedo in the 
visible; $ is the phase function in each regime; B\ is the 
Planck function; and e is the infrared bolometric emissivity. 
T = T(p v q 1 1], e, 9, <f>) is the temperature, which is a func- 
tion of p v ; e; the "beaming parameter," r\; surface planeto- 
graphic coordinates 6 and <p; and the (dimensionless) phase 
integral, q (see below for discussions of r\ and q). 

In practice, the thermal flux depends sensitively on the 
temperature distribution across the surface of the target, and 
uncertainties about that temperature distribution typically 
dominate the uncertainties in the derived albedos and sizes 
(see Fig. 1). Given knowledge of the rotation vector, shape, 
and the distribution of albedo and thermal inertia, it is in 
principle possible to compute the temperature distribution. 
Unsurprisingly, none of these things are known for a typi- 
cal object where we seek to use the radiometric method to 
measure the size and albedo. The usual approach is to use 
a simplified model to compute the temperature distribution 
based on little or no information about the object's rotation 
axis or even rotation period. 

2.1. Standard Thermal Model 

The most commonly employed model for surface tem- 
perature on asteroidal objects is the Standard Thermal 
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Fig. 1. — Thermal models for KBO 38628 Huya (2000 EB173). Spitzer Space Telescope 24 and 70 pm data are shown as circles, with 
vertical error bars within them indicating the measurement uncertainties. Six models are fit to the data, with the resulting model albedos, 
diameters, and beaming parameters summarized in the legend. From top to bottom the models are: 1) Hybrid STM fit to 24 and 70 (im 
data, with 77 as a free parameter (the therm model used here), 2) Hybrid ILM fit to 24 and 70 pm data, 3) Canonical STM (77 = 0.756) 
fit to the 24 pm data, 4) Canonical ILM (77 = 1.0) fit to the 24 pm data, 5) Canonical STM fit to the 70 /im data, 6) Canonical ILM fit to 
the 70 pm data. Note the close agreement of the albedos and sizes for models 1 and 2. Fits to data from one band, using the canonical 
asteroid values for 77, result in much larger uncertainties in the derived parameters, particularly the fits to the 24 pm data. 



Model (STM; cf. Lebofsky and Spencer, 1989, and ref- 
erences therein). The STM assumes a non-rotating (or 
equivalently, zero thermal inertia) spherical object, and 
represents the "hot" end-member to the suite of possible 
temperature distributions. Under STM assumptions, the 
dayside temperature depends only on the angular distance 
from the sub-solar point, 6: T(8) = Tqcos^-^O, and the 
temperature is zero on the night side. The sub-solar point 
temperature T = [(1 - A)S/(r}ea)] 1 l A . Here A = qp v 
is the bolometric albedo, S is the solar constant at the dis- 
tance of the object, and a is the Stefan-Boltzmann constant. 
Even though the STM represents the hottest reasonable dis- 
tribution of surface temperatures for an object in radiative 
equilibrium with sunlight, early studies of the emission 
from asteroids showed that their emission was even hot- 
ter than predicted by the STM (Jones and Morrison, 1974; 
Morrison and Lebofsky, 1979). That led to the introduction 
of the beaming parameter, 77, which allows for localized 
temperature enhancements on the dayside, e.g. in the bot- 
toms of craters or other rough features, and the tendency of 
such warm regions to radiate preferentially in the sunward 
(and, for outer solar system objects, observer-ward) direc- 
tion (i.e. to beam). (Note that while 77 appears analogously 
to the emissivity, e, in the expression for the surface temper- 
ature, 77 does not appear explicitly in the expression for the 
thermal emission, Eq. lb.) Lebofsky et al. (1986) derived a 
value of 77 = 0.756 based on 10 pm observations of Ceres 
and Pallas. We refer to the STM with 77 set to 0.756 as the 
canonical STM. 



2.2. Isothermal Latitude Model (ILM) 

The cold end-member of the suite of plausible temper- 
ature distributions for an object in radiative equilibrium 
with sunlight is the Isothermal Latitude Model (ILM; also 
known as the fast-rotator model). The ILM assumes a 
spherical object illuminated at the equator and rotating very 
quickly (or equivalently, a slowly rotating object with infi- 
nite thermal inertia). The resulting temperature distribution 
depends only on latitude, <f>: T{<j>) = Tqcos x I a 4>, where 
in this case the sub-solar point temperature is given by 
T = [(l-A)5/(7T77ec7)] 1 / 4 . The factor of tt in this expres- 
sion reduces the subsolar point temperature by 33% relative 
to the STM. Because the ILM is characterized by infinite 
thermal inertia, local temperature variations, and therefore 
beaming, are precluded: thus the canonical ILM assumes 
77 = 1. 

2.3. A Hybrid Thermal Model 

Fig. 1 illustrates the problems inherent in using either the 
STM or the ILM to measure the sizes and albedos of KBOs. 
In particular, none of the 4 canonical STM or ILM models 
fit to either the 24 or 70 pm data (4 lower elements in the 
figure legend) match the observed 24:70 pm color. As a re- 
sult, the systematic uncertainties on the albedos and diam- 
eters, depending only on whether the STM or ILM is used, 
are large: pv is uncertain by a factor > 2.5 for the fits to the 
70 pm data, and is uncertain by a factor of > 17 for the fits 
to the 24 pm data. (Note, however, that the relative efficacy 
of these two wavelengths depends on the temperature of the 
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target: if the thermal spectrum peaks near the 24 fim band, 
observations at that wavelength will be considerably more 
effective at constraining the physical properties of the target 
than indicated by this particular example.) However, if the 
beaming parameter, 77, is allowed to be a free parameter of 
the fit (top 2 elements in the figure legend), both the color 
of the thermal emission and its intensity can be matched. 
More importantly, both the STM and ILM give nearly the 
same diameters and albedos with r\ as a free parameter. The 
basic reason for this is that the 24 and 70 [im data provide a 
direct determination of the temperature of the thermal emis- 
sion from the object; equating that color temperature to the 
effective temperature gives a direct estimate of the size of 
the target, independent of the details of an assumed temper- 
ature distribution (and independent of the visual brightness 
as well). 

While the beaming parameter was introduced to model 
enhanced localized dayside temperatures and infrared 
beaming, it can also mimic the effects of other influences on 
the temperature distribution, such as pole orientation (note 
that the emission from a pole-on ILM is indistinguishable 
from the STM), and intermediate rotation rates and ther- 
mal inertias. For example, a rotating body with non-zero 
thermal inertia will have lower dayside temperatures than 
predicted by the STM, but an STM with a value of 77 larger 
than would be supposed based on its surface roughness will 
have a similar color temperature. Likewise, a quickly rotat- 
ing body with a low thermal inertia will have higher dayside 
temperatures than predicted by the ILM, an effect that can 
be mimicked by an ILM with rj < 1. 

Returning to the top two models in the legend of Fig. 1, 
the STM fit results in ry = 1.09, suggesting that the tem- 
perature distribution on the target (the KBO 36828 Huya) 
is cooler than predicted by the canonical STM with rj = 
0.756. Likewise, for the ILM rj = 0.41, suggesting that the 
surface is significantly hotter than would be predicted by 
the canonical ILM with rj = 1. 

2.4. Thermal Model: Application 

In the following we adopt a thermal model in which the 
beaming parameter, along with size and albedo, are free pa- 
rameters which we use to simultaneously fit observed flux 
densities at two thermal wavelengths, and the constraint 
imposed by the visual brightness of the object. Because 
such models have temperature distributions intermediate 
between the canonical STM and ILM, they can be thought 
of as a hybrid between the two. Further, because the sys- 
tematic uncertainties in the model albedos and diameters 
associated with the choice of hybrid STM or hybrid ILM are 
fairly small relative to the uncertainties in the measured flux 
densities and other model assumptions, we simply adopt the 
hybrid STM as our model of choice. (The error bars on py 
and D stemming from the choice of STM or ILM hybrid 
model in Fig. 1 are < 4% and < 2%.) We note that a 
number of studies have employed a similar approach with 
variable rj {e.g. Harris, 1998; Delbo et ah, 2003; Fernandez 



etal, 2003). 

In order to use the STM, we must make some assump- 
tions regarding the nature of the thermal emission and vis- 
ible scattering. We assume a gray emissivity, e = 0.9. 
The infrared phase function, (j) ir — O.Olmag/deg, depends 
only weakly on the emission angle. For our observations, 
emission angles for all but 5 targets (29P, Asbolus, Elatus, 
Thereus and Okyrhoe) were < 5 deg. Because the effects 
are small relative to other uncertainties in the models and 
data, we have neglected the IR phase effect for all of the 
results presented here.. We assume standard scattering be- 
havior for the the objects in the visible, i.e. a scattering 
assymetry parameter, G = 0.15, leading to a phase inte- 
gral q = 0.39 (Bowell et ah, 1989). This assumption also 
allows us to directly relate the geometric albedo py, the 
diameter D, and the absolute visual magnitude, Hy via 
D = 1346 p|/ 2 10- Hv / 5 , where D is in km (Bowell et 
al, 1989; Harris, 1998). By utilizing the absolute visual 
magnitude in this way, the scattering phase function, $„j S 
apparently drops out; however, if the actual scattering be- 
havior differs from the assumption above, our albedos and 
diameters will still be affected because the scattering be- 
havior determines the value of q. We note, also, the results 
of Romanishin and Tegler (2005), who found that absolute 
magnitudes available through the IAU Minor Planet Center 
and through the Horizons service at the Jet Propulsion Lab- 
oratory have are biased downward (brighter) by 0.3 magni- 
tudes. The Hy values shown in Table 1 are culled from the 
photometric literature, and should be fairly reliable. 

For low albedo objects, the albedos we derive depend 
only weakly on the assumed value of q, while for high- 
albedo objects the value of q exerts a strong influence (see 
expressions for Tq in Sections 2.1 and 2.2). For the exam- 
ple of 38628 Huya (Fig. 1), changing to q = 0.8 makes 
only a < 1% difference in the albedo. However, if we use 
q = 0.39 to model the data for the 4 largest objects in the 
sample, 90377 Sedna, 136199 Eris, 136108 (2003 EL 6 i), 
and 136472 (2005 FYg), we obtain geometric albedos that 
exceed a value of 2. While not (necessarily) unphysical, 
such high values for the geometric albedo are unprece- 
dented. Pluto's phase integral q = 0.8, so for these 4 objects 
(only) we adopt that value instead. 

2.5. Thermophysical Models 

More sophisticated extensions to the STM and ILM in- 
clude the effects of surface roughness and (non-zero, non- 
infinite) thermal inertia (Spencer, 1990), and viewing ge- 
ometries that depart significantly from zero phase (Harris, 
1998). However, for the purpose of determining KBO albe- 
dos and diameters from their thermal emission, the hybrid 
STM gives results and uncertainties that are very similar to 
those obtained through application of such thermophyscal 
models (e.g. Stansberry et al, 2006). Because the hybrid 
STM is much simpler, and it produces results comparable to 
thermophysical models, we employ only the hybrid STM. 
(We note that thermophysical models are of significant in- 
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terest for objects where the pole orientation and rotational 
period of the target are known, because such models can 
then constrain the thermal inertia, which is of interest in its 
own right). 

3. SPITZER OBSERVATIONS 

Roughly 310 hours of time on the Spitzer have been al- 
located to attempts to detect thermal emission from KBOs 
and Centaurs, with the goal of measuring their albedos and 
diameters. Spitzer has a complement of three instruments, 
providing imaging capability from 3.6 - 160 fira, and low- 
resolution spectroscopy from 5 - 100 ^m (Werner et al, 
2004). The long-wavelength imager, MIPS (Multiband 
Imaging Photometer for Spitzer, Rieke et al, 2004), has 24, 
70 and 160 ^m channels. Because of the placement of these 
channels, and the sensitivity of the arrays (which are at least 
10 times more sensitive than previous far-infrared satellites 
such as IRAS and ISO), MIPS is well-suited to studying the 
thermal emission from KBOs. 

3.1. The Sample 

Spitzer has targeted over 70 KBOs and Centaurs with 
MIPS. About 2/3 of the observations have been succesful 
at detecting the thermal emission of the target, although 
in some of those cases the detections have a low signal- 
to-noise ratio (SNR). Here we describe observations of 47 
KBOs and Centaurs made during the first 3 years of the mis- 
sion, focusing on observations of the intrinsically brightest 
objects (i.e. those with the smallest absolute magnitudes, 
Hy), and of the Centaur objects. Table 1 summarizes the 
orbital and photometric properties of the sample. 

The distribution of the objects in terms of dynamical 
class is also given, in two forms. The second to last col- 
umn, labeled "TNO?", indicates whether the orbital semi- 
major axis is larger than Neptune's. By that measure, 31 
of the objects are trans-Neptunian Objects (TNOs), and 17 
are what might classically be called Centaur objects; that 
classification is nominally in agreement with the classifica- 
tion scheme proposed in the Gladman et al. chapter, al- 
though they classify Okyrhoe and Echeclus as Jupiter fam- 
ily comets, rather than Centaurs. Another classification 
scheme has been proposed by Elliot et al. (2005; see also 
Dones et al. chapter) as a part of the Deep Ecliptic Sur- 
vey (DES) study, and the target classification thereunder 
appears as the last column in the table. According to the 
DES classification, 21 of the targets in the Spitzer sample 
are Centaurs. 

Thus, about 30^-0% of the sample we discuss here are 
Centaurs, and the rest KBOs. Among the KBOs, only 4 
objects are Classical, while 12 are in mean-motion reso- 
nances with Neptune, 9 are in the scattered disk, and one 
(90377 Sedna) is in the extended scattered disk: Classical 
objects are under-represented. Because Classicals do not 
approach the Sun as closely as the Resonant and Scattered 
Disk objects, and because they have somewhat fainter ab- 
solute magnitudes, the Classicals are at the edge of Spitzer 



capabilities. One Spitzer program has specifically targeted 
15 of the Classicals, but data analysis is ongoing. 

The visible photometric properties of the sample are di- 
verse, and generally span the range of observed variation 
except in terms of the absolute magnitudes, which for the 
KBOs are generally Hy < 7. The spectral properties 
of KBOs and Centaurs are reviewed in the chapters by 
Barucci et al, Tegler et al, and Doressoundiram et al: 
here we summarize those characteristics as regards our sam- 
ple. The visible colors, given in Table 1 as the spectral 
slope (measured relative to V), cover the range from neu- 
tral to very red (Pholus). Visible absorption features have 
been reported in the 0.6-0.75 /jm region for 47932 (2000 
GN m , 38628 Huya, and (2003 AZ 84 ) (Lazzarin et al, 
2003; de Bergh et al, 2004; Fornasier et al. 2004). Sev- 
eral of the targets exhibit near-IR spectral features, with wa- 
ter and methane ices being the dominant absorbers iden- 
tified. Water ice detections have been made for 10199 
Chariklo, 83982 (2002 GO g ), 47171 (1999 TC 36 ), 47932 
(2000 GN m ), 90482 Orcus, 50000 Quaoar, and 136108 
(2003 EL 6 i). 55638 (2002 VE 95 ) exhibits methanol ab- 
sorption, as does 5145 Pholus, along with its strong water 
ice absorption. Methane ice is clearly present on 136199 
Eris (Brown et al 2005), and 136472 (2005 FY 9 ) (Licandro 
et al, 2006a; Tegler et al, 2007; Brown et al, 2007), and 
136199 Eris may also have N 2 ice (Licandro et al. 2006b). 
Two of the objects exhibit surface heterogeneity: 31824 
Elatus (Bauer et al, 2003) and 32532 Thereus (Barucci et 
al, 2002; Merlin et al, 2005). 

3.2. The Observations 

Most of the targets in the sample presented here were 
observed in both the 24 and 70 fim channels of MIPS. In 
a few cases, when the target was predicted to be too faint 
to observe in the second channel, only one channel was 
used. Integration times vary significantly, ranging from 200 
- 4000 sec. As Spitzer observations of KBOs and Cen- 
taurs proceeded, it became clear that they were significantly 
harder to detect than had been predicted prior to the launch. 
The difficulty was due to a combination of worse than pre- 
dicted sensitivity for the 70 array (by a factor of about 
2), and the fact that KBOs are colder and smaller than as- 
sumed. As these realities made themselves evident, later 
observing programs implemented more aggressive observ- 
ing strategies, and have generally been more successful than 
the early observations. 

In some cases the same target was observed more than 
once. These observations fall into three categories: re- 
peat observations seeking to achieve higher sensitivity 
(e.g. 15875 (1996 TP 66 ) and 28978 Ixion), multiple vis- 
its to characterize lightcurves (20000 Varuna and 47932 
(2000 GN171) are the only cases, and neither observation 
produced a measurable lightcurve), and multiple visits to 
allow for the subtraction of background objects (so-called 
"shadow observations"). The basic idea of a shadow ob- 
servation is to observe the target, wait for it to move out of 
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Fig. 2. — Processing of the Spitzer 24 pm data for 90482 Orcus. The left panel shows the typical quality of image available through 
the data pipeline, with scattered light and dark latent artifacts still present. The center panel shows the improvements that can be made 
by correcting the aforementioned artifacts, and reflects the quality of the data we analyzed for targets that were imaged only once. The 
right panel shows further improvement due to the subtraction of a shadow observation, and reflects the quality of data we analyzed for 
targets that were imaged two or more times. 



the way, then re-observe the field. By subtracting the two 
images, the emission from stationary sources is removed. 
Fig. 2 illustrates the shadow method, as well as some of 
the extra processing we apply to the Spitzer 24 fim data to 
improve its quality. 

3.3. Photometry 

Flux densities were measured using aperture photome- 
try, as described in Cruikshank et al. (2005) and Stansberry 
et al. (2006). The apertures used encompassed the core of 
the PSF, out to about the first Airy minimum (their angular 
radii were 10" and 15" at 24 and 70 /mi). Small apertures 
were used to maximize the signal-to-noise ratio (SNR) of 
the measurements. Sky measurements were made in the 
standard way, with an annulus surrounding the object aper- 
ture, and also by placing multiple circular apertures in the 
region around the target when the presence of background 
sources or cirrus structure dictated. The photometry was 
aperture corrected as described in Engelbracht et al. (2007) 
and Gordon et al. (2007). Finally, we apply color correc- 
tions to our measurements as described in Stansberry et al. 
(2007), resulting in monochromatic flux densities at the ef- 
fective wavelengths of the 24 and 70 /mi filters (23.68 and 
71.42 /im, respectively). The MIPS calibration is defined in 
such a way that the color corrections for stellar spectra are 
unity. Even though our targets are much colder (typically at 
or below 80 K), the ~ 20% passbands of the MIPS filters 
result in color corrections that are typically less than 10%. 

Uncertainties on the absolute calibration of MIPS are 
4% and 5% at 24 and 70 /im, respectively (Engelbracht 
et al, 2007; Gordon et al, 2007). In our photometry of 
KBOs and Centaurs we adopt systematic uncertainties of 
5% and 10%, to account for the absolute calibration un- 
certainty and additional uncertainties that may be present, 
e.g., in our aperture and color corrections. At 70 /im our 



adopted systematic uncertainty includes significant margin 
to account for degraded repeatability for faint sources. Ad- 
ditional uncertainty comes from the finite SNR of the detec- 
tions themselves, which is estimated from the statistics of 
values falling in the sky annulus and/or sky apertures. We 
root-sum-square combine the systematic uncertainty with 
the measurement uncertainty determined from the images 
to estimate the final error bars on our measurements, and 
use those total uncertainties in estimating the physical pa- 
rameters we report. The SNR values we tabulate below re- 
flect the errors estimated from the images, and so provide 
an estimate of the statistical significance of each detection. 

4. SPITZER RESULTS 

Our flux density measurements, and the albedos and di- 
ameters we derive from them, are given in Tables 2 and 3. 
Table 2 gives our results for those objects observed in both 
the 24 and 70 /im channel. When only an upper limit on 
the flux density was achieved, the results in Table 2 bound 
the albedo and diameter of the target. Table 3 gives the re- 
sults for those objects observed at only one wavelength, and 
gives a second interpretation of the data for those objects in 
Table 2 that were only detected at one wavelength. 

In both Tables 2 and 3 we give the color corrected 
flux density of each target, the SNR of the detections, and 
the temperature we used to perform the color corrections. 
Where we did not detect the source, we give the 3cr up- 
per limit on the flux density, and the SNR column is blank. 
When an object was not observed in one of the bands (Ta- 
ble 3 only), the flux and SNR columns are blank. In both 
Table 2 and 3, albedos (pv), diameters (D), and beaming 
parameters (77) follow the fluxes and temperatures. 
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4.1. Two- Wavelength Results 

As discussed earlier and demonstrated in Fig. 1, the 
model-dependent uncertainties in the albedo and diameter 
we derive for targets detected at both 24 and 70 pm are 
much smaller than those uncertainties for objects detected 
in only one of those bands, and in particular are usually very 
much smaller than for objects detected only at 24 /im. For 
this reason, we focus first on the targets we either detected 
at both wavelengths, or for which we have constraints on 
the flux density at both. We use these results to inform our 
models for targets with single-band detections and limits. 

We apply the hybrid STM to the observed flux densities 
as follows. For targets detected in both bands (Table 2), 
we fit the observed flux densities and the la error bars, de- 
riving albedo and diameter values and la uncertainties on 
them. For those objects with an upper limit in one band 
and a detection in the other, we fit the detection and the 
the upper limit in order to quantitatively interpret the con- 
straints the limit implies for the albedo and diameter. For 
this second class of observation, we also perform a single- 
wavelength analysis (see Table 3) in order to derive inde- 
pendent constraints on these properties. While the results 
given in Table 2 include values of the beaming parameter, 
rj, those values only reflect the departures of the measured 
emission from the assumptions of the STM; had we chosen 
to model the data with the ILM, the fitted values for rj would 
be entirely different (even though py and D would be very 
similar). Results from observations made at very similar 
epochs are averaged. An exception to that rule is the two 
observations of 38628 Huya. Those data were analyzed in- 
dependently to provide a check on the repeatability of our 
overall data analysis and modeling methods for a "bright" 
KBO, and show agreement at the 4% level for py, and at 
the 2% level for D. 

The average behavior of the targets is of particular inter- 
est for interpreting single-wavelelength observations, where 
we have no independent means for constraining r\. Restrict- 
ing our attention to those targets detected at SNR> 5 at both 
24 and 70 /mi, and excluding the highest and lowest albedo 
object from each class, we find that for outer solar system 
objects the average beaming parameter is rj = 1.2 ± 0.35. 
We re-examine the average properties of the sample later. 

4.2. Single- Wavelength Results 

Because we are primarily interested in the albedos and 
sizes of our targets, we fit our single-wavelength observa- 
tions with the STM, setting the beaming parameter to the 
average value determined above: we term this model the 
"KBO-tuned" STM. We also apply the canonical STM and 
ILM (i.e. with r\ — 0.756 and 1.0, respectively) to the 
single-wavelength data, to interpret the data in the context 
of these end-member models and assess the resulting uncer- 
tainties in model parameters. 

Table 3 gives the results for the single-wavelength sam- 
ple, including those objects in Table 2 with a detection at 
one wavelength and an upper-limit at the other. Where a 



model violates a flux limit, the corresponding albedo and di- 
ameter entries appear as a "?". The albedos and diameters 
we derive using the average beaming parameter from the 
two- wavelength sample are in the columns labeled "KBO- 
Tuned STM"; the range of albedos and diameters resulting 
from application of the canonical STM and ILM are labeled 
"STM " and "ILM ". Note that the flux densities for ob- 
jects in both Table 2 and 3 are sometimes slightly differ- 
ent, because in Table 3 the color correction is based on the 
blackbody temperature at the object's distance, rather than 
on the 24:70 /im color temperature. 

4.3. Spitzer Albedos and Diameters 

The results presented above include low SNR detections, 
non-detections, and multiple results for some targets. In the 
top portion of Table 4 we present results for the 39 targets 
that were detected at SNR > 5 at one or both wavelengths. 
The results for targets that were visited multiple times are 
averaged unless one observation shows some indication of 
a problem. Targets with an upper limit in either band ap- 
pear in both Tables 2 and 3; in the top portion of Table 4 
we give values that are representative of all of the earlier 
models. The top portion of the table contains 39 objects, 26 
detected at both 24 and 70 /im, 9 at 24 /im only, and 4 at 
70 /im only. 17 of the objects have orbital semimajor axes 
inside Neptune, and 21 exterior to Neptune's orbit. Where 
other albedo and diameter determinations exist, the table 
summarizes the result, the basis of the determination, and 
the publication. 

4.4. Other Constraints on py and D 

The albedos and sizes of about 20 TNOs several Cen- 
taurs have been determined by other groups using various 
methods; the lower portion of Table 4 presents those results 
not given in the top portion of the table, and the constraints 
that can be derived from Spitzer data, when those exist (al- 
though the SNR for all 5 cases is low, and for 90377 Sedna 
only a 70 /im limit is available). 

In general our results and those of other groups agree at 
the < 2a level (e.g. 10199 Chariklo, 26308 (1998 SMi 65 ), 
47171 1999 TC 36 , 55565 2002 AWi 97 , 136199 Eris, 
136108 (2003 EL 6 i)). In a few cases there are discrepan- 
cies. For example, our results for 20000 Varuna are incon- 
sistent with the millimeter results of Jewitt et al. (2001) and 
Lellouch et al. (2002), which suggest a significantly larger 
size and lower albedo. While our detection at 70 /im nomi- 
nally satisfied the 5a threshold for Table 4, the background 
showed significant structure and the SNR of the detection in 
the individual visits was actually quite low. Combined with 
the fact that we were not able to directly fit the beaming pa- 
rameter, we are inclined to favor the submillimeter results 
for this object over those from Spitzer. While there is some 
tendency for the Spitzer diameters to be smaller and albe- 
dos higher, there is generally good agreement between our 
Spitzer results and those from other groups and methods. 
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5. ALBEDO STATISTICS and CORRELATIONS 

The Kuiper Belt is full of complexity, in terms of the 
dynamical history and the spectral character of its inhab- 
itants. It is natural to look for relationships between the 
albedos of KBOs and their orbital and other physical param- 
eters. Fig. 3 shows the Spitzer albedos for detections with 
SNR> 5 (top portion of Table 4) as a function of orbital 
semimajor axis, a, perihelion distance, q & , object diame- 
ter, D, and visible spectral slope, S. Because of their sig- 
nificant intrinsic interest, the data for 136108 (2003 ELgi) 
and 90377 Sedna are also plotted. Immediately appar- 
ent in all of these plots is the marked distinction between 
the largest objects (136199 Eris, 136108 (2003 EL 6 i), and 
136472 (2005 FY 9 )) and the rest of the objects. 90377 
Sedna probably also belong to this class, although our data 
only place a lower bound on its albedo. 136199 Eris and 
136472 (2005 FY 9 ) both have abundant CH 4 ice on their 
surfaces, and so are expected to have very high albedos. 
90377 Sedna's near-IR spectrum also shows evidence for 
CH4 and N 2 ices (Barucci et al, 2005; Emery et al, 2007), 
and Schaller et al (2007) show that those ices should not 
be depleted by Jean's escape: it seems likely 90377 Sedna's 
albedo is quite high. The surface of 136108 (2003 EL 6 i) is 
dominated by water ice absorptions, with no evidence for 
CH4 or N 2 , yet also has a very high albedo. Charon, which 
has a similar spectrum, has py — 37%, but some Satur- 
nian satellites (notably Enceladus and Tethys) have albedos 
> 80% (Morrison et al, 1986). 

The dichotomy between 136199 Eris, 136108 (2003 EL 6 i), 
136472 (2005 FY 9 ), Pluto (and probably 90377 Sedna) and 
the rest of the KBOs and the Centaurs, in terms of their 
albedos and spectral characteristics, suggest that they are 
members of a unique physical class within the Kuiper Belt 
population (see chapter by Brown et al). We will refer to 
these objects as "planetoids" in the following, and generally 
exclude them from our discussion of albedo statistics and 
correlations because of their obviously unique character. 

5.1. Albedo Statistics 

Table 5 summarizes the statistics of the Spitzer-dehved 
albedos, and the correlations between albedo and other pa- 
rameters. Because there is no clearly preferred way to dif- 
ferentiate Centaurs from KBOs, we give results for two def- 
initions: a < 30.066 AU (which we term the MFC Defini- 
tion, referring to the Minor Planet Center classification (see 
the Gladmann et al chapter), and the DES Definition (re- 
ferring to the Deep Ecliptic Survey classification (Elliot et 
al. 2005; Dones et al chapter). 

Typical geometric albedos for all of the KBOs and Cen- 
taurs are in the range 6.9%-8.0%, depending on whether the 
mean or median is used, with a dispersion of about 4.1%. 
Regardless of which Centaur classification one chooses, it 
appears that Centaurs may have slightly lower albedos than 
KBOs, although the differences are not statistically signif- 
icant relative the the dispersion of the albedos within the 
classes. The Kuiper variant of the Kolmogorov-Smirnov 



(K-S) test gives no evidence that the albedos of the KBOs 
and Centaurs are drawn from different parent populations, 
regardless of whether the MPC or DES definition of Cen- 
taur is used. Typical values for the beaming parameter 
(exluding results based on an assumed beaming parameter) 
are in the range 1.1 - 1.20, with a dispersion of about 0.4. 
This is in good agreement with the value of 1.2 ± 0.35 we 
adopted for the "KBO Tuned STM" used to construct Ta- 
ble 3. There does not appear to be any significant difference 
in the beaming parameter between KBOs and Centaurs. 

5.2. Albedo Correlations 

Because our errors are non-symmetric and probably non- 
gaussian, we apply the Spearman rank-correlation test to as- 
sess the significance of any correlations between albedo and 
other parameters. Table 5 gives the correlation coefficients 
(p), and their significance (x) in standard deviations from 
the non-correlated case. The albedos for the 4 planetoids 
mentioned above are not included in these calculations. 

Fig. 3a and 3b show pv as a function of the orbital prop- 
erties a and q® . There is an upward trend of py vs a, with 
the objects at a < 20 AU clustering at py ~ 5%, while at 
larger distances there is significantly more scatter in py. As 
shown in Table 5, the correlation between py and a for the 
entire sample is significant at the \ = 2.7a level (99.4% 
likelihood). It appears that most of the correlation is due 
to the Centaurs, but the significance of the Centaur correla- 
tion depends considerably on which definition of Centaur is 
used. (Note that because the number of objects in the KBO 
and Centaur subsamples is about half that of the full sam- 
ple, the significance of the correlations for the subsamples 
is typically lower than that for the full sample.) Because 
the significance of the py vs a correlation is below 3a, it 
is tentatitive. Another reason to treat the correlation with 
some skepticism is that it could reflect biases in the param- 
eter space for KBO discoveries: low-albedo objects will be 
harder to detect at visible wavelengths, and the difficulty 
increases significantly with distance. Because our sample 
is drawn from optically discovered objects, one might ex- 
pect a trend such as seen in Fig. 3a even if there is no real 
correlation between p v and a. 

Fig. 3b reveals a similar correlation between p v and 
<7©, and Table 5 suggests that in this case the correlation 
is significant at the \ = 3.5er level (99.95% likelihood). 
This correlation holds up fairly well for both Centaurs and 
KBOs, regardless of which classification is used. It is pos- 
sible that this correlation could also be due to the discov- 
ery bias mentioned above. However, if it reflects an ac- 
tual relationship between py and q, there may be a fairly 
simple explanation. Objects closer to the Sun will tend 
to experience higher temperatures, depleting their surfaces 
of volatile molecules (which typically have high visible re- 
flectances). Likewise, UV-photolosis and Solar wind radi- 
olysis will also proceed more quickly closer to the Sun, and 
could darken those surfaces (although radiolysis by cosmic 
rays probably dominates beyond about 45 AU (see Cooper 
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Fig. 3. — Geometric albedo plotted vs. a) orbital semimajor axis, b) orbital perihelion distance, c) object diameter, and d) the slope of 
the object's visible spectrum (i.e. color). Open circles are for Centaur objects (a < 30.066 AU); filled circles are for TNOs. In panel 
a) the point for 90377 Sedna has been plotted at a — 150 AU rather than at its true semimajor axis of 489 AU. In panel b) the point for 
90377 Sedna has been plotted at q = 50 AU, rather than at its true perihelion distance of 78 AU. 
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etal. chapter). 

Fig. 3c and 3d show pv as a function of intrinsic proper- 
ties of the objects: the diameter, D, and the visible spectral 
slope (color), S. Fig. 3c shows an apparent correlation be- 
tween pv and D, and particularly so for the KBOs. This 
correlation is apparently confirmed in Table 5, where for 
the MPC classification the pv vs. D correlation is signif- 
icant at the x = 3.4cr level (99.9% likelihood). However, 
for the DES classification the significance is only x = 2.8cr 
(99.5% likelihood), so the correlation is not robust against 
small changes in which objects are considered as KBOs. 
Including the planetoids in the correlation calculation in- 
creases the significance of the correlation to well above 3a, 
but doing so results in a (probably) false impression that the 
albedos of all KBOs are well correlated with diameter. At 
this time it is difficult to conclude that any such correlation 
exists at a statistically significant level. 

Fig. 3d shows an apparent correlation between pv and 

5, particularly for the Centaurs. Table 5 shows that this cor- 
relation is the second most significant for a subclass, with 
2-6 < x < 2.9 (depending on the classification chosen), 
second only to the py vs. D correlation for KBOs. Here, 
the Kuiper variant K-S test does indicate a high likelihood 
(99.95%) that the albedos of red KBOs and Centaurs (with 
S > 0.2) are drawn from a different parent population than 
the gray ones, a similar result to that found based on the 
Centaur colors alone (see Tegler et al. chapter). A nat- 
ural assumption might be that the color diversity of KBOs 
and Centaurs results from mixing between icy (bright, spec- 
trally neutral) and organic (dark, red) components. How- 
ever, this correlation suggests that red objects systemati- 
cally have higher albedos than the gray ones. On the basis 
of spectral mixing models between spectrally neutral dark 
materials (such as charcoal) and red material (represented 
by Titan tholin), Grundy and Stansberry (2003) suggested 
that just such a correlation between red color and higher 
albedo might exist. Why the Centaurs might embody this 
effect more strongly than the KBOs is still a mystery. Inter- 
estingly, the three most spectrally neutral objects defy the 
color-albedo trend, having rather high albedos: there may 
be at least two mechanisms underlying the observed color 
diversity. Those objects are 2060 Chiron, 90482 Orcus and 
(2003 AZ 84 ), and their unique position in the albedo-color 
plane may indicate that they share some unique surface 
character. 

6. FUTURE PROSPECTS 

At present, Spitzer/MlPS provides the most sensitive 
method available for measuring thermal fluxes from typi- 
cal KBOs, but several upcoming observatories and instru- 
ments will provide substantially improved sensitivity. The 
joint ESA/NASA Herschel mission will have at least a fac- 
tor of 2 better sensitivity at 75 microns (compared to MIPS 
70 micron sensitivity), and additionally have a number of 
photometry channels in the range 70-500 microns. Since 
cold KBOs have their thermal emission peaks in the range 



60-100 microns, observations in the Herschel bandpasses 
will map the peak of a KBO SED. Herschel is scheduled 
for launch in late 2008. The Large Millimeter Telescope 
(LMT) in central Mexico will have sufficient sensitivity at 
1 millimeter with the SPEED instrument to detect thermal 
flux from the Rayleigh-Jeans tail of cold KBOs. First light 
for the LMT is expected in 2008. 

Farther in the future, the American-European-Chilean 
Atacama Large Millimeter Array (ALMA) will provide suf- 
ficient sensitivity from 0.35-3 millimeters to detect typical 
KBOs; first light for ALMA might be as soon as 2012. The 
Cornell Caltech Atacama Telescope (CCAT) will operate 
at 200 microns to 1 mm, and its sensitivity at 350 yum will 
surpass that of ALMA; first light could also be in 2012. 
Any of Herschel, ALMA, and CCAT (the case is less con- 
vincing for the LMT) could be used for a large survey of 
many moderate-size (100 km class) KBOs. Such a pro- 
gram would expand the number of KBOs with good ther- 
mal measurements (and therefore radii and albedos) from 
tens to hundreds. 

All of these next-generation capabilities operate at wave- 
lengths either near the emission peak of KBOs, or well out 
on the Rayleigh-Jeans part of their spectra. While albe- 
dos and diameters derived from such observations are less 
model-dependent than those based on single-wavelength 
observations taken shortward of the emission peak, there 
are still significant uncertainties. For example, canonical 
STM and ILM fits to an 850 fim flux density produce albe- 
dos that differ by about 30%; if the KBO-tuned STM is 
used (including its uncertainty on 77), that uncertainty is cut 
almost in half. If the validity of the KBO-tuned STM is 
born-out by further Spitzer observations of KBOs, it can be 
used to significantly refine the albedos and diameters de- 
rived from sub-millimeter KBO detections. 

7. SUMMARY 

Efforts to characterize the physical properties of KBOs 
and Centaurs with Spitzer are beginning to pay off. Con- 
siderable improvements have been made in the first three 
years of the mission in terms of predicting the necessary 
integration times, developing aggressive and successful ob- 
serving strategies, and data processing. We present our 24 
and 70 /im observations for 47 targets (3 1 with orbital semi- 
major axes larger than that of Neptune, 16 inside Neptune's 
orbit), and apply a modified version of the Standard Ther- 
mal Model to derive albedos and diameters for them. 39 of 
the targets were detected at signal-to-noise ratios > 5 at one 
or both wavelengths. We use that sample to look for rela- 
tionships between albedo and the orbital and physical pa- 
rameters of the objects. The most marked such relationship 
is the distinct discontinuity in albedo at a diameter of about 
1000 km, with objects larger than that having albedos in 
excess of 60%, and those smaller than that having albedos 
below about 25%. We suggest that these large, very high 
albedo objects (90377 Sedna, 136108 (2003 EL 6 i), 136199 
Erisand 136472 (2005 FYg)) constitute a distinct class in 
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terms of their physical properties. 

The data suggest possible correlations of albedo with or- 
bital distance, and with size and color, but the statistical 
significance of the correlations is marginal. Two correla- 
tions, those of albedo with perihelion distance (for KBOs 
and Centaurs) and with diameter (for KBOs), are nominally 
significant at more than the 3cr level. Perhaps the most inter- 
esting trend (albeit significant at only about the 2.8a level) 
is for distinctly red Centaurs to have higher albedos than 
those that are more gray, contrary to what might intuitively 
be expected. 

Prospects for improving on and expanding these results 
are relatively good. Spitzer will be operational into 2009, 
and more KBO observations will probably be approved. 
New ground- and space-based observatories will also con- 
tribute significantly, and at wavelengths that are comple- 
mentary to those used here. In particular, submillimeter- 
millimeter studies of KBOs should be relatively easy with 
facitilties such as ALMA, CCAT and LMT The Herschel 
mission should also be very productive at far-IR to submil- 
limeter wavelengths. 
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a Small body number, provisional designation, and proper name for the target sample. 
b Orbital semimajor axis (a), eccentricity (e) and inclination (i). 

c Absolute Visual Magnitude (H v ), and spectral slope and uncertainty (S and erg, in % per 100 nm relative to V band), 
from the photometric literature. 

d Orbital semimajor axis > that of Neptune (30.066 AU). 

e Deep Ecliptic Survey dynamical classification (Elliot et at, 2005): CENTR = Centaur, CLSCL = Classical, RESNT = 
Resonant, SCTNR = Scattered Near, SCEXT = Scattered Extended. 
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Two-Band Thermal Model Results 
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9041920 


18.816 


18.350 


3.528 


78.0 


29.7 


6.1 


76.0 


44 +1.27 
• 3 - 44 -0.82 


206.7;^.! 


1 fiq-°- 30 

j-'^+o.so 


55565 


(2002 AW197) 


9043712 


47.131 


46.701 


0.155 


7.7 


15.0 


6.7 


51.9 


_ +4.42 
1L -< '-3. 00 




1 26~ ' 20 
i - ZD +0.22 


55576 


Amycus (2002 GBi ) 


17766144 


15.589 


15.155 


6.367 


86.1 


13.6 


5.8 


99.9° 


1 7 Qfi +7 ' 77 


7fi S — 12 ' 5 


64"° 18 
u -°^+0.19 


55636 


(2002 TX 300 ) 


10676992 


40.979 


40.729 


<0.065 




<11.1 




48.4 


> 17.26+f 3 f 


< 641.2^|||| 


2 16 -0 ' 78 
z - lo +0.95 


55637 


(2002 UX 25 ) 


10677504 


42.368 


42.413 


0.486 


15.0 


23.0 


5.3 


57.2 


11.0U_ 3 og 


681.2 +115 ' 6 


1 04"° 18 
i - u4 +0.18 


60558 


Echeclus (2000 EC 98 ) 


8808960 


14.141 


13.736 


4.901 


84.7 


15.5 


5.0 


94.0 


, oq+189 
•J-o-J-l.OS 




, 9 c-0.32 
1 - zd +0.33 


65489 


Ceto (2003 FXi2 8 ) 


17763840 


27.991 


27.674 


1.463 


71.5 


14.6 


12.2 


73.6 


7 fJ7+ 1 ' 38 

'■"'-1.10 


229-7;^. 6 2 


0.86;°;?? 



TABLE 2 — Continued 



Number 1 Name (Designation)* AORKEY b Rq c A c F 24 d SNR 2 4 d F 7 o d SNR 70 d T 24 :70 e Pv l (%) -D 



73480 


(2002 PN 34 ) 


17762816 


14.608 


14.153 


10.368 


99.0 


31.0 


12.6 


95.3 


„ __i_n 83 
4 - 25 -0.65 


- - „ . — in 2 
119 - 5 +10.3 


. 1n -0 15 


83982 


Crantor (2002 GOg) 


9044224 


14.319 


13.824 


2.276 


58.6 


<8.7 




>89.8 


> 8-60±3.36 


< 66.7+™;? 


-1 a a — 57 


84522 


(2002TC302) 


13126912 


47.741 


47.654 


0.054 


6.5 


18.0 


3.1 


44.8 


3. 08IJ 24 


1 1 a<~ a — 325.0 
114.5-4+337.4 


Q qq— 0.54 

Z - d6 +0.53 


84922 


(2003 VS 2 ) 


10680064 


36.430 


36.527 


0.304 


6.0 


25.7 


3.5 


52.8 


5.84+4-^ 

— 2 24 


79c; 9 — 187.6 
+199.0 


2.007n 

+0.54 


90482 


Orcus (2004 DW) 


13000448 


47.677 


47.442 


0.329 


32.4 


26.6 


12.5 


53.1 


19.72+1 4 7 ° 


946.3^_ 7 ^ ^ 


i n«-° 09 


90568 


(2004 GV 9 ) 


13000960 


38.992 


39.007 


0.166 


18.2 


17.5 


9.2 


51.4 






1 94- 020 

1,ayi +0.20 


119951 


(2002 KX14) 


10678016 


39.585 


39.197 


<0.109 




<11.7 




51.2 


> 8.09±| |f 


< 561.6+J?J;g 


1 01-0.66 

1 ' al +0.84 


120061 


(2003 COi) 


17764864 


10.927 


10.917 


21.722 


99.0 


33.4 


11.3 


114.7 


5 74+1.49 


76.9+11 


u ' yl +0.20 


136108 


(2003EL 6 i) g 


13803008 


51.244 


50.920 


<0.022 




7.8 


5.3 


<44.6 


? h 


? h 


? h 


136199 


Eris (2003 UB 3 i 3 )s 


15909632 


96.907 


96.411 


<0.014 




2.7 


4.0 


<40.1 


? h 


? h 


? h 


136472 


(2005 FY 9 ) s 


13803776 


51.884 


51.879 


0.296 


21.1 


14.6 


9.4 


54.8 


? h 


? h 


9h 




(2002 MS 4 ) 


10678528 


47.402 


47.488 


0.391 


20.5 


20.0 


5.1 


56.6 


°- 4i -2.26 


79fi o - 1 22 - 9 


o.88; °:i« 




(2003 AZ 84 ) 


10679040 


45.669 


45.218 


0.291 


12.4 


17.8 


6.7 


55.2 


12 S2+ 4 ' 31 


685.8 +98 8 





a Small body number, provisional designation, and proper name for the target sample. 
b Unique key identifying the data in the Spitzer data archive. 
c Target distance from the Sun and Spitzer, in AU. 

d Color-corrected flux densities (mJy) at 23.68 fim and 71.42 /im. Upper limits are 3<r. SNR is signal to noise ratio in the images (see text). 

°The temperature of the blackbody spectrum used to compute the color correction. In most cases this is the 24:70 £tm color temperature, but for the 4 denoted targets, the subsolar 
blackbody temperature was lower than the color temperature, and we used that instead. 

f The visible geometric albedo (p v , percentage), diameter (D, km) and beaming parameter (77) from hybrid STM fits. Fits to upper limits provide a quantitative interpretation of the 
constraints they place on py and D. 

^Results for 136199 Eris, 136108 (2003 EL 6 i) and 136472 (2005 FY 9 ) assumed a phase integral of 0.8, typical of Pluto. 

h No STM with plausible albedo and beaming parameter can simultaneously fit the 24 and 70 (im data. For 136472, models with two albedo terrains can fit the data, and give D ~ 1500 
km. 



Table 3 

Single-Band Thermal Model Results 























KBO- 


-Tuned STM 


Pv" (%) 




Number a 


Name (Designation) 41 


AORKEY a 




A a 




SATi? 2 4 a 


-F70 a 


SAffl 70 a 


b 

T.ss 


Pv (%f 


D c 


STM ILMo 


STM ILMo 


5145 


Pholus (1992 AD) 


9040896 


18.614 


18.152 


3.119 


66.0 


<19.6 




91.4 


8.wt 6 ? W 


138.6;, 340 


17.07 - ? 


95.8 - ? 


5145 


Pholus (1992 AD) 


12661760 


19.827 


19.768 


0.987 


18.8 


<10.0 




88.6 




QQ 4 — 2 3.2 
-1-25 


32.74-? 


69.2 - ? 


7066 


Nessus(1993 HA 2 ) 


9033984 


19.501 


19.219 


0.440 


12.4 






89.3 


6.53^2 46 




14.02-1.44 


40.8-127.4 


10370 


Hylonome (1995 DW 2 ) 


9038080 


19.963 


19.824 


0.530 


14.9 


<10.0 




88.3 






13.28-1.32 


47.9- 152.0 


10370 


Hylonome(1995DW 2 ) 


12659968 


20.333 


20.390 


0.451 


16.0 






87.5 


6.33 +5 2 f 2 


69.3-^ | 


13.80- 1.34 


46.9-150.5 


15875 


(1996 TP 66 ) 


8805632 


26.491 


26.250 


0.720 


17.9 


<17.3 




76.6 




191 S" 54 ' 9 


12.54-? 


123.1 -? 


15875 


(1996 TP 66 ) 


12659456 


26.629 


26.113 


0.437 


14.6 


<6.8 




76.4 


8.2lt 7 / 61 


152.2; 42 '' 6 


19.37-? 


99.1 -? 


20000 


Varuna(2000 WR 106 ) 


9045760 


43.209 


42.830 


<0.094 




10.9 


4.9 


60.0 


? ' 


? 


? - 8.09 


? -744.1 


20000 


Varuna(2000WRioe) 


9031680 


43.261 


43.030 






10.0 


5.6 


60.0 


17.77^ 


502.o; 6 6 4 9 ;= 


26.34-8.68 


412.3-718.2 


28978 


Ixion(2001 KX 76 ) 


12659712 


42.510 


42.058 


0.303 


7.9 


<18.3 




60.5 


25.81 d 


446.3 d 


32.28-? 


399.1 -? 


31824 


Elatus(1999 UG 5 ) 


9043200 


10.333 


9.998 


5.990 


69.8 


<12.3 




122.7 


6.417I?- 52 


41.3-?' 1 


11.41-? 


31.0-? 


31824 


Elatus(1999 UG 5 ) 


12661248 


11.125 


10.826 


8.596 


99.0 


<8.9 




118.3 


? 


? 


? - ? 


? - ? 


35671 


(1998 SN1135) 


9040384 


37.967 


37.542 






14.7 


6.3 


64.0 




458.2;^ 1 1 


6.42-2.17 


376.4-648.1 


42355 


Typhon (2002 CR 46 ) 


9029120 


17.581 


17.675 






31.4 


8.6 


94.1 


5 09+ 1 ' 24 


1T3.8+J|;° 


6.81-3.13 


150.3-221.7 


52975 


Cyllarus (1998 TF 35 ) 


9046528 


21.277 


21.001 


0.274 


8.7 






85.5 


11 46+ 8 - 96 


61-9^iI.I 


24.43 - 2.42 


42.4-134.9 


63252 


(2001 BL 41 ) 


9032960 


9.856 


9.850 


4.864 


95.6 






125.7 


3.90^-^ 


34-2; 6 6 ; 5 7 


6.93-1.34 


25.6-58.3 


83982 


Crantor(2002 G0 9 ) 


9044224 


14.319 


13.824 


2.310 


58.6 


<8.6 




104.2 


11.18+™ 9 


58.5; ? 12 ' 7 


21.28-? 


42.4 - ? 


90377 


Sedna(2003 VBi 2 )° f 


8804608 


89.527 


89.291 






<2.4 




41.7 


> 20.9li^ 


< 1268.8; 2 g 2 'J 


32.93-8.17 


1010.9-2029.0 


136108 


(2003 EL 6 i)° 


13803008 


51.244 


50.920 


<0.025 




7.7 


5.3 


55.1 


84.11±»;« 


ii5i.o;| 9 ; 9 


96.41 -59.12 


1075.1 - 1372.9 


136199 


Eris (2003 UB 3 i 3 ) e 


15909632 


96.907 


96.411 


<0.014 




2.7 


4.0 


40.1 


68.91™ 


2657.0; 20 «;« 


84.90-39.17 


2393.7-3523.9 


136472 


(2005FY 9 )° 


13803776 


51.884 


51.879 


s 




14.6 


9.4 


54.8 


70 on +10.30 
78.20_ 8 55 


1502.9; 9 "; 2 


91.63-52.55 


1388.3 - 1833.3 


136472 


(2005 FY 9 ) c 


13803776 


51.884 


51.879 


0.296 


21.1 


s 




54.8 


oo.ya_ 12 25 


2215.2; 399 ; 2 


59.34-6.27 


1725.3-5307.0 



a The first 9 columns are identical to those in Table 2. Flux densities that are blank indicate no data exist. 

b The subsolar temperature of a blackbody at the distance of the target. Color corrections are made using a black body spectrum with this temperature. 

c The range of visible geometric albedos (given as a percentage) and diameter (in km) derived from fitting the KBO-tuned STM (i.e. r) = 1.2 ± 0.35), and the canonical STM and ILM. "?" indicates the the model emission 
violates a flux limit. 

d Only the KBO-tuned STM using r\ — 0.85 did not violate the 70 fim flux limit for this observation of Ixion. 

"Results for 90377 Sedna, 136108 (2003 EL 6 i), 136199 Eris and 136472 (2005 FY 9 ) assumed a phase integral of 0.8, typical of Pluto. 
1 Fit to the 70 fim upper limit: lower bound on p v , upper bound on D. 

g Fits to the individual bands for 136472 (2005 FYg) are shown: it is not possible to simultaneously fit both bands with a single thermal model. 



Table 4 
Adopted Physical Properties 



Number a 



Name (Designation) 2 



Physical Properties from Spitzer Data 



^detect" TNO? a Method c 



Other Methods 



29P Schwassmann-Wachmann 1 

2060 Chiron (1977 UB) 



5145 Pholus(1992AD) 

7066 Nessus (1993 HA 2 ) 

8405 Asbolus (1995 GO) 

10199 Chariklo(1997CU 26 ) 

10370 Hylonome (1995 DW 2 ) 

15875 (1996 TP 66 ) 

20000 Varuna (2000 WRi 06 ) 

26308 (1998 SMi 65 ) 

26375 (1999 DE 9 ) 

28978 Ixion (2001 KX 76 ) 

29981 (1999 TDio) 

31824 Elatus (1999 UG 5 ) 

32532 Thereus (2001 PT13) 

35671 (1998 SNi 65 ) 

38628 Huya (2000 EB173) 

42355 Typhon (2002 CR 4 e) 

47171 (1999 TC 36 ) 

47932 (2OOOGN171) 

50000 Quaoar (2002 LM 60 ) 

52872 Okyrhoe (1998 SG35) 

52975 Cyllarus(1998TF 35 ) 

54598 Bienor (2000 QC243) 

55565 (2002 AW197) 

55576 Amycus (2002 GBi ) 

55637 (2002 UX 25 ) 

60558 Echeclus (2000 EC 9g ) 

63252 (2001 BL41) 

65489 Ceto (2003 FXi 28 ) 

73480 (2002 PN34) 

83982 Crantor (2002 GO9) 



4.61 



+5.2 2 
1.90 



- --+1.03 

1 - ol -0.87 



8.0 



+7 
3 



6.5 
5.46 



+5.3 
2.5 
1.27 



- 7 o+0.49 



+5 

-3 

7A± 7 3 
16+i° 



6.33 
6.85 
12 
4.40 
10 
4.28 
4.3; 



1.53 
1.16 
1.58 
1.19 
+ 14 



hi. 42 
-0.96 

hi. 09 
-0.80 
-1.8 

5-04+8:!? 
+1.3 
0.9 
1.53 

-1.17 



5.1 
7.18 



5.68 



+2.54 
1.59 

13.2 



2.49] 
11.5" 1 



19.91 7 
vfO.'si 

-0.55 
-+9 

3 44 +i.27 

77 +4.42 

17 96+ 7 ' 77 

11 50+ 5 ' 09 
n.au_3 og 

o oo+l. 89 
' J - o * > -1.08 

3.9±™ 
7.67t\f 
4-25±g;|| 



37.3 



233.3 



+ 11.3 

14.4 



+ 14.7 



140 



+40 



60 



84.2 



+16 
7.8 



+7.8 



258.6 



+10.3 



70 



160 



+20 
45 



500 



+45 
-100 
+ 100 



279.8^i;? 
461.0T 453 



650 
103.7 
30 
87.8 
460 
532.6 
175 
414.6 



+46.1 

2 20 
+260 
13.5 



+8 

9.4 
+9.5 
80 
+60 

24.4 
+25.1 

20 
+ 17 

38.2 
+38.8 



321.0 



844 
52.1 



+57.4 

190 
+207 
-6.9 
l +6.9 

62~ 18 
DZ +20 

206 7~ 30 1 
ZUD -'+30.1 
7S4 fi~ 108 ' 3 

' D - J +12.5 
681.2 +n5 ' 6 



35 
229.7 
119.5 



60 



+7 
18.2 
18.6 
-10.2 
hl0.3 
13 



0.26 
1.13 



+0.28 
-0.13 
+0.14 



1.3 



+0.4 



1.2 
0.80 

1.23 

1.2 
1.2 
1.2 

1.48 
1.05 
0.8 



0.35 
+0.35 
0.16 



+0.10 

0.35 
+0.35 
0.35 
+0.35 
0.35 
+0.35 



1.64 
1.2 



+0.17 
0.12 
+0.12 
0.2 
+0.2 
-0.31 
+0.32 
0.35 
+0.35 

1 50~ ' 28 
± - ou +0.30 
-0.35 
'+0.35 

1 - uy +0.07 
1 o -0.35 
1 - z +0.35 

1 - L '+0.13 



1.27 



+ 15 



r) 09 — U.IO 

z - c>z +0.46 

1 4" ' 4 
1 -*+0.4 

1 -*°+0.39 
1 r)-0.35 
i - z +0.35 

1 69~ ' 30 
1 - oa +0.30 

1 26~ ' 20 
x - zu +0.22 

64" ' 18 

1 04~ ' 18 
1 - u ^+0.18 
1 oc-0.32 
1 - zd +0.33 

1 q-0.35 
1 - z +0.35 

n «fi-°' 09 
U - 8b +0.10 
1 j-0.16 
1 - 1 +0.16 
1 20 — °' 35 



both 


Cen 


mIR 


13 ±4 


both 


Cen 


mTR 


17 ±2 






ISO 


11 ±2 






mIR 


14 ±5 


24 


Cen 


mTR 


4.4+1.3 






IRS 


7.2 ± 2 


24 


Cen 






both 


Cen 


mIR 


12 + 3 






IRS 


4.3 ± 1.4 


both 


Cen 


mm 


5.5 ± 0.5 






mIR/mm 


7 ± 1 


24 


Cen 






24 


TNO 






70 


TNO 


submm 


6 ± 2 






mm 


7 ± 3 


both 


TWO 

11NU 


mm/bin 


y. 1 ± 4 


both 


TNO 






24 


TNO 


mm 


> 15 


both 


TNO 


IRS 


6.5 


24 


Cen 


IRS 


5.7 ± 2 


both 


Cen 






70 


TNO 






both 


TNO 


mm 


> 8 


70 


TNO 






both 


TNO 


mm 


5 ± 1 






mm/bin 


14 It 


both 


TNO 






both 


TNO 


image 


9±3 


both 


Cen 






24 


Cen 






both 


Cen 






both 


TNO 


mm 


9±2 


both 


Cen 






both 


TNO 






both 


Cen 






24 


Cen 






both 


TNO 






both 


TNO 






24 


Cen 







40 ± 5 CrS3 
144 ± 8 Fe02 
142 ± 10 Gn05 

180 Ca94 
189 ± 26 Da93 
148 ± 25 Cr06 



66 ± 4 Fe02 
95 ± 7 Cr06 
275 AW2 
246 ± 12 Gn05 



1016 ± l56 Je01 ' A(04 

914 ± 156Le02,Al04 

238 ± 55 Ma04 - G !' 05 

< 804 Ai04 

98 Cr06 
36 ± 8 Cr06 



< 540 AW4 
609 ± 70 AW4 

302 ± 7 QMa04,Gy05 

1260 ± 190 Br04 



977 ± 130 Ma02 



TABLE 4 — Continued 



Physical Properties from Spitzer Data Other Methods 



Number a 


Name (Designation) 21 






if- 


\ b 

A <2etect 


TNO? a 


Method 




D a 


yu4sz 


( /TOD/I T~*\X7\ 

wrcus tzuu4 tjw ) 


1Q 79 +3.40 
19.72_2. 76 


O/ifi Q — 72.3 


l-U8 +0 . 10 


both 










90568 


(2004 GV9) 


8.051! .46 


677 - 2 +71.3 


1 n/i— 0.20 
L94 +0.20 


both 


TNO 








120061 


(2003 COi) 


- » . -1- 1 4Q 


76.9 + |:| 


0-91+Sio 


both 


Cen 








136108 


(2003 EL 6 i) 


84.I20 1 


1150.^™ 


70 


TNO 


Lcurve 


65 ±6 


1350 ± 100 Ra05 


136472 


(2005 FY 9 ) 


80+2°; 


1 500 ~ 200 
iouu. +4 q 




both 


TNO 










(2002 MS 4 ) 


8.41+j;™ 


VOfi O — 122.9 
,ZD - Z +123.2 


0-88 + o;j4 


both 


TNO 










(2003 AZ 84 ) 




685.8 +98 ; 8 


1 04 -016 


both 


TNO 








15789 


14993 SO 










TNO 


ISO 


3.5 ± 1.4 


298 + 140 1 " 00 


15874 


(1996TL fifi ) 


3.5+?-° 

— 1.1 


°'°+116 


1 8" - 3 

+0.3 


both 


TNO 


ISO 


> 1.8 


< 958 ™° 


19308 


H996 TfWl 










TNO 




> 3.3 


< 902 A '0i,Gy05 


19521 


Chaos (1998 WH24) 










TNO 


mm 


> 5.8 


< 747 Ai04 > G !/05 


24835 


(1995 SM55) 










TNO 


mm 


> 6.7 


< 704 a;o4,G!/05 


55636 


(2002TX 300 ) 


> 10 


< 800 




limit 


TNO 


mm 


> 19 


< 709 Or04,Gy<)5 


58534 


(1997 CQ29) 










TNO 


bin 


39 ±17 


77 ± 18 Ma04 > No04,Gy05 


66652 


(1999 RZ253) 










TNO 


bin 


29 ± 12 


170 ± 39 N °04,Gy05 


84522 


(2002TC302) 


i+2.9 


1 1 ^n -325 
11JU +33r 


r, n-0.5 

z - J +0.5 


both 


TNO 


mm 


> 5.1 


< 1211 Ai04 . G y° 5 


88611 


(2001 QT297) 








TNO 


bin 


10 ±4 


168 ± 38 Os03 ' G y° 5 


90377 


Sedna (2003 VB i2 ) 


> 16. 


< 1600. 




limit 


TNO 


image 


> 8.5 


< 1800 Br04a 


136199 


Eris (2003 UB313) 


70 + 15 - 


2600.; 2 °° 




70 


TNO 


mm 


60 ±8 


3000 ± 200 Be06 














image 


86 ± 7 


2400 ± 100 Br06 




(1998 WW31) 










TNO 


bin 


6 ±2.6 


152 ± 35 Ve02 ' G « 05 




(2001 QC298) 










TNO 


bin 


2.5 ± 1.1 


244 ± 55Ma04,Gj/05 



NOTE. — Results above the horizontal line have Spitzer detections with SNR > 5; those below the line have SNR < 5, or no Spitzer data. 
a Columns 1-5 and 7 are as defined in Table 2. 

b Wavelengths where the objects were detected at SNR> 5 (above horizontal line), or have lower quality Spitzer data (below line). 

c "Method" by which the diameter was measured. The meanings are: "bin" (binary mass plus density assumption), "image" (HST upper limit), 
"IRS" (Spitzer mid-IR spectra), "ISO" (Infrared Space Observatory), "Lcurve" (lightcurve + rotation dynamics), "mIR" (Groundbased 10-20 ^m), 
"mm" (typically 1.2 mm groundbased data), "submm" (typically 850 /im groundbased data) 

References. — (A102) Altenhoff et al. (2002); (A104) Altenhoff et al. (2004); (Be06) Bertoldi et al. (2006); (Br04) Brown and Trujillo (2004); 
(Br04a) Brown et al. (2004); (Br06) Brown et al. (2006); (Ca94) Campins et al. (1994); (Cr83) Cruikshank and Brown (1983); (Cr06) Cruikshank et 
al. (2006); (Da93) Davies et al. (1993); (Fe02) Fernandez et al. (2002); (Gn05) Groussin et al. (2004); (Gy05) Grundy et al. (2005); (JeOl) Jewitt et 
al. (2001); (LeOl) Lellouch et al. (2002); (Ma02) Margot et al. (2002); (Ma04) Margot et al. (2004); (No04) Noll et al. (2004); (Or04) Ortiz et al. 
(2004); (Os03) Osip et al. (2003); (Ra05) Rabinowitz et al. (2005); (ThOO) Thomas et al. (2000); (Ve02) Veillet et al. (2002) 



Table 5 
Geometric Albedo 



MPC Classification 21 DES Classification 13 







All 


Centaurs 


KBOs 


Centaurs 


KBOs 


Quantity 








Statistics 






Avgerage 




8.01 


6.55 


8.87 


6.30 


9.88 


Median 




6.85 


5.74 


7.67 


5.73 


8.41 


a c 




4.07 


2.68 


4.22 


2.50 


4.23 


#Obj. d 




35 


15 


18 


19 


14 










Correlations 






Parameter 


P e 


x f 


P e x f 


p £ x f 


P e x f 


P e x f 


a 


0.46 


2.74 


0.70 2.78 


0.24 1.03 


0.41 1.81 


0.16 0.60 




0.58 


3.49 


0.58 2.32 


0.65 2.83 


0.43 1.94 


0.53 2.04 


D 


0.45 


2.70 


-0.08 0.32 


0.77 3.37 


-0.07 0.31 


0.72 2.80 


S 


0.40 


2.40 


0.64 2.58 


0.13 0.56 


0.66 2.94 


-0.08 0.32 



a Centaurs classified as objects having orbital semimajor axes < 30.066 AU. 

b Centaurs classified by dynamical simulations (Deep Ecliptic Survey, Elliot et al. (2005)). 

c Standard deviation of the albedo values. 

d Number of Spitzer albedos used (from Table 4). The highest and lowest values were excluded. 
e Spearman rank correlation coefficient between albedo and the parameter in the left column. 
'Significance of the correlation, in standard deviations relative to the null hypothesis. 
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